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Abstract 


The Institute for Energy and Environment (IEE) at the University of Strathclyde has developed various fuel cell (FC) systems for stationary 
and vehicular applications. In particular the author is involved in the development of alkaline fuel cell (AFC) systems. To understand the dynamic 
behaviour of the system’s key element, the alkaline fuel cell stack, a dynamic model was developed allowing the characterisation of the electro- 
chemical parameters. The model is used to forecast the behaviour of the fuel cell stack under various dynamic operating conditions. The so-called 
Nernst potential, which describes the open circuit voltage of the stack, is calculated using thermodynamic theory. Electrochemistry theory has 
been used to model the sources of the electric losses within the FC, such as activation, ohmic and concentration losses. The achievable value of 
this paper is the first publication of a detailed dynamic AFC based on mass balance, thermodynamics and electrochemical theory. The effects of 
the load changes on various fuel cell parameters, such as electrolyte concentration and concentrations of dissolved hydrogen and oxygen were 
covered in this investigation using the author’s model. The model allows a detailed understanding of the dynamic effects within the AFC during 


load change events, which lead to the experienced electric response of the overall FC stack. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The Centre of Economic Renewable Power Delivery 
(CERPD) at the University of Strathclyde has developed various 
fuel cell (FC) systems for stationary and vehicular applications 
over the last 5 years. The aim of the research is the design and 
build of reliable and cost efficient hybrid FC/battery systems, 
which could replace existing conventional technology in the near 
future. A domestic scale combined heat and power (CHP) alka- 
line fuel cell (AFC) system has been developed. The small AFC 
stack (1-3 kW.) is used to satisfy the average load demand of the 
domestic load profile of a stationary application. Whereas the 
load peaks are supplied by a battery system working in parallel 
with the AFC stack. 

However, using an AFC for a highly fluctuating domestic 
load needs precise knowledge of the dynamic behaviour of the 
fuel cell system and consequently of its main component, the 
fuel cell stack. This paper describes a dynamic electrochemi- 
cal model of an alkaline fuel cell stack. The practical electric 
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losses of the fuel cell, such as activation, ohmic and concentra- 
tion losses are considered and the variations of internal fuel cell 
parameters, such as reactant concentrations and pressures, dur- 
ing steady-state operation and dynamic load changes are shown 
and explained. 


2. AFC operation 


An AFC operates by introducing hydrogen and oxygen (in 
case of the modelled Zetek stack provided as air) gases into the 
gas diffusion layers of the anode and cathode, respectively. The 
gases subsequently diffuse into the catalyst layers of the elec- 
trodes, where they partly dissolve into the KOH electrolyte. To 
avoid gases entering the electrolyte layer the potassium hydrox- 
ide solution (KOH) in the electrolyte compartment has a slightly 
higher pressure than the gases in the anode and cathode. The 
normal operating pressure of hydrogen and air for the FC were 
set to 40 mbar above atmospheric pressure, whereas the KOH 
pressure is set to 100-150 mbar. The oxidation reaction occur- 
ring in the catalytic layer of the anode causes the hydrogen to 
decompose and yields electrons to the current collector and the 
remaining protons react with the hydroxy] ions available in the 
KOH solution [1]. 
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Electrode variables 


e Material (Ni, PTFE) 
e Surface area (al, a8) 


e Geometry (T,€) 
e Surface condition 


Electrochemical 

variables 

* Exchange current 
densities (i,) 

¢ Transfer coefficient (a) 

e Limiting current 


density (i,) 


Anode 


Mass transfer variables 


e Mode: (diffusion, convection, 
migration) 

e Flux (N; ) 

+ Diffusitivity (D;) 

e Henry constants (H;) 


Solution variables 


* Bulk concentration of dissolved 
electroactive species (C®,,,,C®,,) 


e Surface concentration of dissolved 
electroactive species (C%,,,,C%5) 


Electrical variables 


e Potential (E) 

e Current (1) 

e Resistance (R) 

* Conductivity (0) 
e Capacitance (C4) 


Operating variables 


¢ Temperature (T) 
+ Pressure (p,.p,) 
¢ Time (t) 


Structural variables 
Cathode 


+ Number of cell (Ngc) 

e Connection of cells 
(parallel, series) 

« Layer thickness (x,) 

e Cell size (Age, 


Fig. 1. Variables affecting the dynamic behaviour of an AFC stack. 


At the cathode, the catalytic reduction of oxygen with water 
of the electrolyte and electrons supplied from the anode, to gen- 
erate hydroxyl ions. The electro-osmotic drag within the cell, the 
electrolyte concentration gradient and the movement of charge 
particles in the electric field, moves the hydroxyl ions towards 
the anode. The electrode reactions as explained are given in the 
following equations: 


Hə +2(OH~) —> 2H20 + 2e7 (1) 
H2O + 0.50) +2e7 > 2(OH7) (2) 


Hence, the overall reaction of the fuel cell can be expressed 
as 


H2 + 0.502 > H20 (3) 


3. The AFC stack model 


The AFC behaviour is complex and is influenced by a large 
number of parameters. The model parameters can be grouped 
into four categories: 


(1) Operating parameters: such as gas inlet pressures and tem- 
perature, reference electrolyte concentration. 

(2) Electrochemical parameters: such as exchange current den- 
sity of the electrodes, transfer coefficients of oxidation and 
reduction processes. 

(3) Material parameters: such as composition of the electrodes 
and the current collectors, porosity of the materials, con- 
ductivity of electrodes. 

(4) Structural parameters: such as layer thicknesses, active area 
sizes, number of cells. 


Fig. 1 shows a summary of the most important fuel cell 
variables, which influence the dynamic stack performance. 


The detailed model developed is based on a combination of 
the variables and parameters shown in Fig. 1, which are related 
through mass balance, mass transfer theory, thermodynamics, 
electrochemical and electrical theories [2]. 

For the purpose of analysis the alkaline fuel cell was split 
into five layers, each of which was modelled separately. These 
layers are: 


anode gas diffusion layer; 
anode catalyst layer; 
electrolyte layer; 

cathode catalyst layer; 
cathode gas diffusion layer. 


Within each layer the appropriate state parameters, e.g. the 
partial pressures of the gases within the diffusion layers, the 
concentrations of the reactant species within the catalyst layer, 
were determined. The different layers are linked through input 
and output parameters. The modelled single cell performance 
was scaled according to the number of cells in the stack. It was 
assumed that all cells in the stack behave in a similar way. The 
physical structure of the AFC and a schematic of the top layer 
of the fuel cell model are shown in Figs. 2 and 3, respectively. 

The model has been developed within the simulation package 
Matlab/Simulink. The block diagram shows the sub-systems of 
the anode and cathode gas diffusion layers, the anode and cath- 
ode catalyst layers and the electrolyte layer in the middle. The 
input parameter of the fuel cell stack model is the load current 
demand oag. The main output parameter is the FC stack voltage 
Estack- 


3.1. Model equations 


3.1.1. Gas diffusion layers 
Each electrode possesses a gas diffusion layer, which is nec- 
essary to provide the process gases uniformly and extract the 
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Fig. 2. Layer structure of the alkaline fuel cell. 


product species. They are designed to prevent the electrolyte 
from weeping into the gas stream. Consequently, the model 
assumes that the hydrophobicity of the electrode will prevent 
any liquid from entering this region. Therefore, only hydrogen 
gas and water vapour will exist in the anode gas diffusion layer, 
and oxygen gas, nitrogen gas and water vapour will exist in the 
cathode gas diffusion layer. 

The diffusion gas transport through porous material is 
described by the one-dimensional Stefan—Maxwell diffusion 


LEVEL 2 


Anode Catalyst Layer 


Air in 
Water out 
L ere rla- rja L ccol >| 
3 4 5 
equation: 
dy; RT 
Be ONN (4) 


g 
j pDij 


The list of symbols is given in Appendix A. To describe the 
mass balance in the gas diffusion layer the continuity equation 
has been used. The general one-dimensional continuity equation 
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Fig. 3. The detailed alkaline fuel cell model. 
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is given by the following equation: 


1 deSPL p; dN; 
= + RP + R° 5 
RT dt dx : : (6) 


Within the gas diffusion layer there occurs no chemical reac- 
tion and no mass transport across a phase boundary (from gas to 
liquid or reverse), hence R? and Rf equals zero. The continuity 
equation becomes: 


1 del’: dN; 


RT dt dx (6) 


The anode contains two species, hydrogen and water, whereas 
the cathode is assumed to contain oxygen, nitrogen and water. 


3.1.2. Catalyst layers 

The mass balances on anode and cathode are defined by the 
one-dimensional continuity equation. In the catalyst layers the 
reaction gases hydrogen and oxygen dissolve into the liquid 
electrolyte. The rate that the hydrogen and oxygen gases dissolve 
in the electrolyte is expressed by the mass transport rate over 
the phase boundary RP . In equilibrium the transport of species 
i across the phase boundary, i.e. the amount of gases dissolved 
in the liquid electrolyte, can be described as follows 


dnp Hipi — Ci 
R= i— aD} ( 2e ‘) (7) 


dx 


The amount of oxygen and hydrogen that react in the catalyst 
layer is expressed with the electrochemical reaction rate RẸ. 
The reaction interface area is not equal to the contact area of the 
CL and electrolyte layer; it is rather defined by the area where 
the electrolyte is in contact with the active material within the 
catalyst layer. This reaction area is expressed through a! the 
specific area of catalyst—electrolyte interface. The reaction rate 
can be represented with Faraday’s law: 

dN¢& siahiloe 


RE = = 8 
i dx ziF 8) 


The gas transport is modelled using the Stefan—Maxwell 
equation for gas diffusion in porous material (Eq. (5)) and the 
mass transport of the liquid species was described by the one- 
dimensional Nernst—Plank equation: 


dC; d® 
N} = —D} i zu; FC;—- — Cu (9) 
dx dx 


The reactants on anode and cathode diffuse, in gaseous form, 
through the ‘dry’ part of the catalyst layer, dissolve in the liquid 
electrolyte and react as stated in Eqs. (1) and (2) on the interface 
area within the catalyst layer. The electrochemical potentials of 
anode and cathode have been described by the Nernst equation. 
Included in the model is the used current—overpotential equation 
which takes into consideration the electric losses under load 
conditions (including the transport losses): 

Si 
RT Ci 
ziF n(&) | my 


EN = EN _ — In 


Gai Gy 
‘loc __ loc II i aaFn/ RT I j —acFn/RT 
L = ıl e = e 
0 ref ref 
Ci j C $ 


(1) 


In the model Eq. (11) was expressed in the form of the Tafel 
equation. Thus, the anode and cathode overpotentials are respec- 


tively: 
qi 
2.303 RT Ci 
Te l loc =| -loc ae l A 
= IT (one s e I(EN T) 
(12) 
2.303RT [ag i Co” 
Ne = ne logie" — logig.. + log II E 
i j 
(13) 


Additionally the double-layer capacitance Cy affects the elec- 
tric response of the fuel cell. The double-layer capacitance of 
the anode and cathode Ca,a, were separately modelled with a 
feed-back loop using the following equation: 
dij, i 1 

dt 7 Ca,a,c 7 Ra,cCa,a,c 


(ni = nic) (14) 


The model calculates the ohmic losses within the catalyst 
layer. Two sources of ohmic losses have been considered: 


e The potential drop within the electrolyte solution defined by 
the solution conductivity. 

e The voltage drop within the electrodes defined by the ohmic 
resistance of the solid electrode materials and the current 
collector. 


The solution potential drop is described by 


1X; 
Ads =" (15) 
Oo 


The voltage drop within the catalyst solid material is 
described by Ohm’s law: 


iAFC 


(16) 


Nohmic,a,c = 
Oa,cXi 
As aconsequence of the electric and electrochemical consid- 
erations the overall anode and cathode voltage was calculated 
from 


Fac = EN, + Nac + a,c + Nohmic,a,c (17) 


3.1.3. Electrolyte layer 

Within the electrolyte layer there are no chemical reactions. 
Further, no gassing or dissolution of reaction species occurs, i.e. 
there is no mass transfer over the phase boundary. Five species 
exist within the layer: dissolved hydrogen, dissolved oxygen, 
water, hydroxyl anions and potassium cations. All species occur 
in the liquid phase. As can be seen from the simulation results, it 
is important to consider the change in electrolyte concentration 
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across the fuel cell, as the electrolyte concentration affects most 
model parameters. 

In the model the electrolyte concentration in the middle of 
the electrolyte layer Ce_2,5 was assumed to be constant and 
equals the reference concentration Ce rer. The electrolyte con- 
centrations at point 2 (Ce_2) and point 3 (Ce_3) were calculated 
using the continuity equation. No chemical reaction and trans- 
port over a phase boundary occurs in the electrolyte layer, hence 
the continuity equation consists only of Nernst—Plank transport 
term and has the following form: 


adc,  — dNe 
de dx 
d ( Dko ppo (dCe/dx)+zg+ Fags Co(db/dx)+(v™*Ce/dx)) 
a dx 


(18) 


Additionally the solution potential ® was calculated in the 
electrolyte layer and is defined by 


1X; 
A@ = — 
oO 


3.1.4. Overall cell and stack potentials 
The anode and cathode potentials Ea, Ee follow from the 
above considerations: 


Ea = EN — Ø, — a — NiR.a (19) 


Ee = EN — Øo — ne — nine (20) 


The overall single cell potential Ece follows from the differ- 
ence of cathode and anode potential: 


Ecell = Ec = Ea (21) 


The model calculates the stack voltage Estack by multiplying 
the simulated voltage of a single cell Ecen with the numbers of 
cells connected in series (Ngeriesmodules ): 


(22) 


Estack = EceliMseries/ modules 


The input load current Joaq Signal was divided by the number 
of cells connected in parallel in the stack parallel to calculate the 
current demanded of a single cell Ken: 


I load 


Feen = (23) 


"parallel 
To take recognition of the ohmic losses caused by the inter- 
connections the overall stack voltage is reduced by an ohmic 
voltage loss Minter: 


(24) 


Estack,real = Estack — Minter 
4. Simulation results 
Using the operating conditions listed in Table 1 the distri- 


bution of the main AFC parameters across a single cell were 
simulated and compared to test and literature data [3-5]. 


Table 1 

Operating conditions of Zetek AFC stack 

Parameter Symbol Value 
Anode and cathode inlet gas pressures Pas Pe 1.053 atm 
Reference electrolyte concentration Ce,ref 6.6M 
Fuel cell temperature Tre 70°C 


4.1. Steady-state simulations 


Fig. 4 shows the ideal voltage—current characteristic of a two- 
module Zetek stack compared with simulation results. The test 
data was supplied by the manufacturer. The data does not include 
information about the activation and concentration limitations, 
only the open circuit voltage of the stack is given (11.3 V). 

It is clear from Fig. 4 that the simulation accurately describes 
the stack polarisation. The correlation factor Reni was calculated 
to be 96.7% for the given data points. The overall stack resis- 
tance Rgtack in the linear region (25-110 A) was calculated to 
be 25 mQ. The limiting current for the stack determined by the 
model is 140 A (130 mA cm~?). The open circuit voltage is cal- 
culated to be 11.3 V. The nominal operating voltage of the stack 
is 8 V (0.67 V per cell) which leads to a nominal stack current of 
100 A. The generated electric power of the fuel cell stack relative 
to the stack current is also shown in Fig. 5. At the nominal operat- 
ing current of 100 A the stack generates 800 We. After reaching 
a current of 125 A the concentration losses increase exponen- 
tially, due to the lack of reaction oxygen (as later shown). At a 
stack current of 140 A the voltage and consequently the power 
collapses. This is the point where the concentration of oxygen 
on the reaction interface equals zero. 

Additionally the polarisation curve of a single cell of the 
Zetek AFC stack was investigated. The characteristics for anode, 
cathode and the overall single cell are shown in Fig. 6. The Zetek 
stack simulated consist of 24 cells connected in parallel (4 cells) 
and series (12 cells) fashion. Hence, the current of a single cell 
was calculated to be a quarter of the fuel cell stack current. The 
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Fig. 4. Zetek AFC stack voltage—current characteristic, comparison between 
simulation results and manufacturer data. 
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Fig. 5. Zetek AFC stack power—current characteristic, comparison between 
simulation results and manufacturer data. 


single cell open circuit voltage was calculated to be 0.94 V, which 
confirms the manufacturer data. The overall losses of the Zetek 
cell are principally determined by the cathode overpotential with 
a simulated limiting current density of 130 mA cm7? (35 A for 
a single cell, 140 A for the stack). 

Fig. 7 shows the electrolyte distribution throughout a single 
cell relative to the applied current. As expected the simulation 
shows that under open circuit conditions the electrolyte concen- 
tration is uniform throughout the three layers with a reference 
concentration of 6.6 M. When the load is applied the KOH con- 
centration on the anode side decreases and Ce on the cathode 
side increases due to the chemical reactions on the electrodes. It 
can be seen that the main concentration drop occurs within the 
electrolyte layer, whereas within the anode and cathode catalyst 
layers the electrolyte concentration varies only slightly. 

This effect can be explained by the fact that the electrolyte 
layer is 1000 times the thickness of the catalyst layers (elec- 
trolyte layer 1 mm, catalyst layers 0.001 mm). The simulation 
shows that it is important to take recognition of the variation in 


Voltage [V] 


oe wee cea pasnccnonpecenaacepaccescce 


fi 
0 5 10 15 20 25 30 35 40 
Current [A] 


Fig. 6. Voltage—current characteristics of anode, cathode and single cell. 
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Fig. 7. Electrolyte concentration within the Zetek alkaline fuel cell as a function 
of current. 


Ce. Thus, the electrolyte concentration has an important influ- 
ence on a variety of fuel cell properties as previously highlighted. 

The concentration profile of the dissolved hydrogen and oxy- 
gen within the catalyst layers and the electrolyte layer relative 
to the current are shown in Figs. 8 and 9. It is clear from 
Fig. 8 that even at the limiting current of 35 A for a single 
cell, sufficient amount of hydrogen can be provided to the reac- 
tion area. The hydrogen concentration drops to a minimum of 
0.9 x 1077 mol cm? in the catalyst layer. It has been assumed 
that hydrogen travelling through the electrolyte to the cathode 
will immediately react on the interface. Hence, the hydrogen 
concentration on the electrolyte—cathode catalyst layer interface 
equals zero. 

Fig. 9 illustrates that at higher currents the oxygen concen- 
tration is substantially reduced. At the limiting current 7 35 A 
(i; =0.13 A cm7?) the concentration of dissolved oxygen in the 
catalyst drops to zero and consequently any further load increase 
will lead to the collapse of the fuel cell voltage. 


= 
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Fig. 8. Dissolved hydrogen concentration profile within the alkaline fuel cell as 
a function of current. 
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Fig. 9. Dissolved oxygen concentration profile within the alkaline fuel cell as a 
function of current. 


In Figs. 10 and 11 the anode and cathode overpotentials (acti- 
vation and concentration losses) are shown, respectively. The 
positive sign of the anode overpotential results from the fact that 
the Nernst potential of the anode is positive (0.3 V) at the given 
operating conditions, whereas the cathode Nernst potential is 
negative (—0.8 V). The higher overpotential on the cathode is 
the result of its oxygen mass transport limitations and lower 
exchange current density. 

As expected (see also Fig. 6) the cathode overpotential is 
much larger than the anode overpotential and consequently the 
cathode overpotential is the performance-dominating factor. It 
is shown that the total anode overpotential is very small, with 
its maximum at 35 A of 0.027 V. Due to the preferable hydro- 
gen oxidation process on the anode the activation overpotential 
at low current (until 10 V) is calculated with the model to be 
zero. This shows that at limiting current (35 A) there is still 
enough dissolved hydrogen in the catalyst layer to keep the oxi- 
dation process going. The overpotential on the cathode reaches 
a maximum of —0.35 V at the limiting current. At low current 
(0-5 A) the activation polarisation is dominant, whereas at a 
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40 
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Fig. 10. Anode overpotential profile as a function of current. 
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Fig. 11. Cathode overpotential profile as a function of current. 


current over 30 A the concentration polarisation causes an expo- 
nential increase of the losses. The figure also illustrates, that the 
model predicts a collapse of the fuel cell voltage for a current 
higher than 35 A through a dramatic increase in the cathode 
overpotential. 


4.2. Dynamic load switching simulations 


Figs. 12 and 13 show the voltage and current responses of the 
fuel cell stack from open circuit to an applied resistive load of 
0.558 Q, respectively. Both figures compare the test results with 
the simulated dynamic behaviour of the stack. 

The open circuit voltage of the stack at the beginning of the 
test was measured to be 10.7 V. After 0.5 s the resistive 0.558 Q 
load was applied. Immediately following the application of the 
load the voltage drops due to the internal resistance of the stack. 
The exponential decay which follows is caused by a combination 
of mass transport dynamics and the partial discharge of the elec- 
tric double layer. The final steady-state voltage level of the stack 
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Fig. 12. Voltage—time characteristic, comparison between test data and simula- 
tion results for load switching of a 0.558 Q resistive load. 
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Fig. 13. Current-time characteristic, comparison between test data and simula- 
tion results for load switching of a 0.558 resistive load. 


is reached after approximately 1 s and is 6.8 V. The test started at 
open circuit. Directly after the load is applied the current rises to 
19.8 A due to the discharge of some stored electric energy within 
the double-layer capacitance structure. Within approximately 
1s after switching the stack current decreases exponentially to 
a steady-state current value of 13.2 A. The model accurately 
predicts the transient behaviour of the 0.558 © resistive load 
switching test with calculated correlation factors Reni of 98.5% 
for the voltage transient simulation and an Reni value of 98.2% 
for the current transient simulation. 

To investigate the effects of the load switch on the fuel cell 
internal behaviour some parameters were investigated in more 
detailed. The electrolyte concentration has an influence on a 
variety of fuel cell parameters and properties. Hence, was inves- 
tigated to determine the changes across the fuel cell during the 
load switching event. Fig. 14 shows the electrolyte concentra- 
tion across the anode CL, the electrolyte layer and the cathode 
CL. At the beginning of the test the fuel cell operates at open 
circuit and the electrolyte concentration is constant throughout 
the liquid layers at 6.6 M. After the load was applied the elec- 
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Fig. 14. Electrolyte concentration profile for load switching of a 0.558 Q resis- 
tive load. 
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Fig. 15. Dissolved hydrogen concentration profile for load switching of a 
0.558 Q resistive load. 


trolyte concentration changed relative to the fuel cell current. On 
the anode the concentration decreases due to the production of 
water, whereas on the cathode the concentration increases due 
to the increase of hydroxy] ions. 

Fig. 14 shows that the change in electrolyte concentration 
for the simulated switching test of 0.558 Q is minimal. On the 
anode and cathode the maximum change is calculated to be 
4 x 1077 mol cm? (4 x 1074 M). However, a comparison with 
Fig. 7 shows that the simulated electrolyte concentrations have 
not reached their steady-state values for the applied current of 
13.2A. 

The steady-state change of the electrolyte concentration is 
approximately 10 times (4 x 10~° mol cm?) larger than the value 
of Fig. 14. The reason of this difference is caused by the slow 
mass transport within the electrolyte layer. 

In Figs. 15 and 16 the concentration profiles of dissolved 
hydrogen and dissolved oxygen in the anode and cathode catalyst 
layers during load switching are shown, respectively. 

After the load is applied at 0.5 s the concentration drops from 
1.093 x 1077 to 1.068 x 1077 mol cm~? and after a further 1 s 
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Fig. 16. Dissolved oxygen concentration profile for load switching of a 0.558 Q 
resistive load. 
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Fig. 17. Anode overpotential of a single cell for load switching of a 0.558 Q 
resistive load. 


adjust to its new steady-state value of 1.076 x 1077 molcm7?. 
The initial concentration drop is caused by the higher reaction 
rate of oxygen during transient times and is a direct result of 
the current spike shown in Fig. 13. The oxygen concentration 
shows a similar behaviour as Cy,. However, due to the smaller 
amount of oxygen used in the reaction the absolute oxygen con- 
centration drop is smaller (1.1 x 107° mol cm~?) compared with 
the Cy, drop (1.7 x 107° molcm~3). Also it must be noticed 
that the amount of dissolved oxygen (1.36 x 1078 mol cm7?) is 
approximately eight times smaller than the amount of dissolved 
hydrogen (1.093 x 1077 mol cm?) on open circuit. The reason 
is as explained before is the lower solubility of oxygen in the 
electrolyte. 

Figs. 17 and 18 show the anode and cathode overpotentials 
caused by the activation and concentration losses, respectively. 
Both overpotentials are zero at the beginning of the load switch- 
ing simulation, as the fuel cell is operating at open circuit. After 
the load is applied the overpotentials rise immediately mainly 
due to the experienced concentration limitations. It is shown in 
Fig. 17 that the simulation predicts a relative small value for the 
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Fig. 18. Cathode overpotential of a single cell for load switching of a 0.558 Q 
resistive load. 


anode overpotential of 0.042 V after the load switching, whereas 
the overpotential of the cathode increases sharply to —0.26 V 
for a single cell (see Fig. 18). The anode and cathode overpo- 
tentials of the stack were simulated to be 0.126 and —0.78 V, 
respectively, for the test stack. 


5. Conclusions 


The dynamic model forecasts the steady-state and dynamic 
behaviour of the alkaline fuel cell stack. The effects of the 
load changes on various fuel cell parameters, such as elec- 
trolyte concentration and concentrations of dissolved hydrogen 
and oxygen were covered in this investigation. The author’s 
model considers the main loss effects of an AFC; the activa- 
tion, ohmic and concentration losses. The dynamics behaviour 
of the fuel cell process was modelled through dynamic mass 
balances of the reactant species and through the consideration 
of the double-layer capacitance. A comparison of simulation 
results with manufacturer and test data has shown the accuracy 
of the model. 

The model will be part for a larger hybrid fuel cell/battery 
system model and is also the basis for the development of an 
AFC condition-monitoring unit. 


Appendix A. List of symbols 


aÈ specific area of gas—electrolyte interface (cm? cm~+) 

a specific area of catalyst-electrolyte interface 
(cm? cm7?) 

AFC area of the cell (cm?) 


Cac double-layer capacitance of anode and cathode (F) 

Ci concentration of species i in electrolyte solution 
(mol cm~?) 

co reference concentration of specie i (mol cm~?) 

D! diffusion coefficient of species i in liquid phase 
(cm? s7!) 

Diy gas diffusion coefficient of species i in j (cm? s7!) 

Eac anode and cathode potential (V) 

EN, anode and cathode Nernst potential (V) 

EY Nernst or open circuit potential at standard conditions 
(V) 

F Faraday constant (96 485 C mol`!) 

H;i Henry constant for species i (mol cm~? bar” !) 

i current density (A cm~?) 

pee anodic and cathodic local current density (A cm~7) 

iv anodic and cathodic local exchange current density 
(A cm~?) 

Ni molar flux of species i (mol cm~? s7!) 

N! molar flux of species i in liquid phase (mol cm~? s7!) 

p total pressure (bar) 

Pi partial pressure of species i (mol cm~?) 

qij reaction order of species i, j 

R universal molar gas constant (8.31451 J mol`! K7!) 

Rac anode and cathode ohmic resistance ({2) 


RE chemical reaction rate of species i (mol cm~? s7!) 
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mass transport rate over phase boundary of species i 
(mol cm~? s7!) 

stoichiometric coefficient of species i 

temperature (K) 

effective mobility of species i (mol cm? J~! s7!) 
thickness of layer i 

vapour mole fraction of species i 

number of electrons involved in reaction 


Greek letters 


Na,c 
L 
Na,c 
i-1 
Nac 
NiR,a,c 


anodic and cathodic transfer coefficient for reaction 
thickness of electrolyte film (cm) 

porosity of the GDL 

anode and cathode overpotential (V) 

anode and cathode overpotentials at present time (V) 
anode and cathode overpotentials at last time step (V) 
anode and cathode ohmic losses (V) 


Nohmic,a,c anodic and cathodic voltage drop (V) 


ye volume average velocity (cm s7!) 

Oac effective anodic and cathodic conductivity (S cm7!) 
Pac anode or cathode solution potential (V) 
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